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We present a method for calculating any (nested) harmonic sum to arbitrary accuracy for all 
complex values of the argument. The method utilizes the relation between harmonic sums and 
(derivatives of) Hurwitz zeta functions, which allows a harmonic sum to be calculated as an ex- 
pansion valid for large values of its argument. A program for implementing this method is also 
provided. 



Q^". I. INTRODUCTION 

o: 

Probably the most efficient method for evaluating cross sections involving parton distribution functions and/or 
fragmentation functions involves the inverse Mellin transform of the Mellin transformed cross section, because the 
H . DGLAP equation can be solved analytically in Mellin space up to the desired order in the strong coupling as [H. In 
' a minimal subtraction scheme, this evolution is determined by perturbatively calculable, mass independent, timelike 
. and/or spacelike splitting functions respectively, which depend on (nested) harmonic sums [2| in Mellin space [1,0] j 
as do massless coefficient functions Q and other quantities in perturbation theory. Harmonic sums are defined for 
positive integer values of the Mellin space variable N as 



X 



rv. c (m \- («gn(fci))"^ (sgn(fc2))"- ^ (sgn(fc3))"^ ... 

O ^' ni = l '^1 122 = 1 "-2 n3 = l "-3 

D ' 

^ However, in the inverse Mellin transform, N is complex. The numerical evaluation of harmonic sums in this case is 
I— the purpose of this letter, and is intended for important calculations such as the cross section calculations of global 

' fits of parton distribution functions (see e.g. Ref. [1]) and fragmentation functions (see e.g. Ref. @). 

^ , The complex analysis of harmonic sums has been understood for some time M • Algebraic relations among harmonic 

QQ sums, namely their quasi-shuffle algebra [1], worked out to weight 6 in Ref. @, and structural relations [lO|, reduces 

■r^lj- \ the number of harmonic sums that needs to be calculated, and can be used to relate them through polynomials to 
SUN), which fixes their single poles at non-positive integer values of N The analytic continuation is then unique 

I [i^l- The continued form of harmonic sums that the method of this letter relies on is the large |iV|, or asymptotic, 

. expansion 

o 

o : 5'=..'=..'^3....(iv) = EE^£t..3....(i-^+W^. (2) 



(We choose to expand in (IniV -I- 7^;) instead of In since this simplifies explicit expressions somewhat.) The coefR- 

;J] ' cients ^^''''^ ^. consist of rational numbers and special numbers such as ln(2), C„, tt" etc. They may also depend 
on 77, where = 1 if the harmonic sum is continued from even integer values of N and ry = — 1 if continued from 
odd values (see Ref. [H, [13] for detailed studies of 77) . To achieve the expansion in eq. ^ , we extend the method 
introduced in the appendix of Ref. [Tsj for the calculation of S* = 5'_2,i. There this harmonic sum was analytically 
continued by extending the sum over ni in eq. ^ to infinity. This approach was also used in Ref. [ll]. In Ref. [IB], the 
harmonic sum was then calculated as an expansion in inverse powers of N. There are 2 problems with this approach 
that prevent it from being applied to other harmonic sums. Firstly, an intermediate step involved the identification of 
a non convergent series with a finite result, which cannot work for all harmonic sums. Secondly, the analytic contin- 
uation only works when fci ^ 1, since otherwise the harmonic sum is logarithmically singular as A^ — > 00. These two 
problems are solved in the present letter, respectively by identifying certain sums with the Hurwitz zeta functions, 
whose expansion for large complex arguments is well known, and reordering the sum when ki — I such that this 
identification also works in this case. We note that a method for obtaining the large |A^] expansion for harmonic sums 
was also given in Ref. 16] using the isomorphism of harmonic sums to shuffle algebra. 

Our approach provides an alternative to that of Ref. 0, [l3] ; where approximations were given for the harmonic 
sums involved in various perturbative quantities such as the 3-loop anomalous dimension. In these references, the 
harmonic sums are determined semi- analytically by calculating the analytic Mellin transforms of the relevant x-space 
functions which are partly calculated by polynomial approximations using the very accurate minimax method, which 
leads to an accuracy of 10"^" to 10~^^. (In fact, only some harmonic sums were determined this way, since the 
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remaining can then be determined by using the various relations between harmonic sums discussed above, some of 
which are presented in Ref. Here we are concerned with deriving a method for the numerical evaluation of any 
harmonic sum anywhere in complex Mellin space to any desired accuracy. 

The rest of the letter is organized as follows. In section [TTl the basic analytic continuation that allows harmonic 
sums to be calculated for complex arguments is discussed. The improvement to this result that leads to fast, accurate 
calculations via an expansion in the inverse of the complex argument is discussed in section IIIII Section IIVI discusses 
the additional manipulations required when the first index of the harmonic sum is equal to 1 , but the method itself is 
essentially the same. The additional calculation required when the absolute value of the argument is not sufficiently 
large for the expansion to work is given in section |Vl Finally, in section IVTl we summarize our approach in full. The 
appendix contains explicit results for harmonic sums whose weight (the sum of the absolute values of the indices) is 
less than or equal to 3. 



II. ANALYTIC CONTINUATION 



Our approach is most conveniently derived by writing eq. ([T]) in the form 

<? (AT\ (sgn(fci))"i 

ni = l 

Equation ^ can be analytically continued to complex values of N as follows: The sum over ni — l,2,...,iV is 
replaced with a sum over ni ~ 1,2,..., oo, and then the additional contribution due to this replacement is canceled 
by subtracting a second term which is identical except with the sum over ni = N + 1, N + 2, . . . , oo. Then, in this 
second term, N is subtracted from ni so that the sum is over ni = 1, 2, . . . , oo. The result is 

Consequently, for any complex N, all harmonic sums can be determined from the 5*^, being the harmonic sums whose 
depth (the number of indices) is equal to 1. 

Although eq. (|4|), which was also derived in Ref. [H, [l^, is essentially our starting point, as it stands it is 
insufficient for the numerical evaluation of harmonic sums because when ki ^ 1 the sum over ni may converge slowly, 
as mentioned in Ref. [17|, while for fci = 1 the series for = oo in eq. ([3]) does not converge. As we will see, the 
improvement in this letter is the development of a method for obtaining the large \N\ expansion of any harmonic 
sum, which converges fast for sufficiently large |A^|. 

III. CALCULATION FOR LARGE |iV| WHEN fci / 1 

In this section, we give the general procedure for obtaining the large |A^| expansion of any harmonic sum, eq. ([2]), 
from the expansion of harmonic sums of lower depth and weight. This means that all harmonic sums' expansions can 
be obtained once the expansion for depth 1 harmonic sums are known. In fact, it suffices to know only the expansion 
for Si{N) [2], 

Si{N)^i;{N + l)+-fE, (5) 
where ip{N) and its derivatives can be evaluated for large |A^| using 

^(iV + l)^ln(iV) + ^-f:^, (6) 

k=l 

since the expansion for the remaining harmonic sums of depth 1, which are also presented in Ref. can be determined 
from eq. (O using eq. We find 

Sk{N) = -^Sk-i{N) + {k- l)C(fc) for k > 2, (7) 
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so that 



Using the result 



^fc(^) = TT^^'/'^'^'H^ + 1) + C(fc) for k>2. 
{k - 1)1 



n=l 



2k 



Sk 



N 



Sk 



N - 1 



for fc > 1 



(8) 



(9) 



that follows from eq. ((T]), one obtains from eq. ^ the S-k for fc > 1 

S-i{N) ={-lff3iN + l)-\n2, 



S-k{N) = 



(fc-1) 



2k 



for fc > 2 



(10) 



where P{N) = (1/2) [V'((A^ + l)/2) - V(^/2)]. 

We now show how eq. ([4]) can be expanded in 1/N once the expansion in eq. ([2]) for Sk2.k3,...{^) is known. We 
assume fci ^ 1, and will deal with the case fci = 1 in section [TVl The expansion for Ski,k2,k3,---i-^) may then be 
calculated by inserting eq. ([2]) with N ^ ni + N into eq. (|4]) and making the identifications 



y ' 



C(x,iV+l) 



E 



ln™(ni+iV)- 



(m + iV)^ 



(±i)"i 



1 

2x 



iV + 1 1 



2 

m oo 



7V + 1 



(11) 



dx ) E 



(±l)"i 



[ni+NY 

where the Hurwitz zeta function C(a;, + 1) is set equal to its large \N\ expansion 

N 1 ^ B2k {2k + x-2y. 1 



C{x,N + l) = N-"" 



1 ^ 

fe=l 



(2fc)! (x-l)! iV2fe- 



(12) 



with i?2fc the Bernoulli numbers. Note that the second result in eq. pT|) follows from the first. 
Equation (|12p can be derived from the result 



ax,N + i) = 



^(^-i)(iV + 1) 
(-l)^(a;- 1)! ' 



(13) 



which follows from the first result in eq. (fTT|) together with eq. ((8]) and the analytic continuation in eq. ([4]) applied to 
Sk for fc > 2 (see also Ref. (3), and the result that follows from eq. ([5]): 



^(^-i)(7V + l) = (-1)^ 



(a; -2)! (a; - 1)! 



N 



2A^^ 



fc=i 



2fe 



(2fc + 2;-2)! 
(2fc)!7V2fe+^-i 



(14) 



Equation (|12p is continued to all real values of x > 1, noting that ({x, + 1) must be real in this case for real values 
ofiV>-l. 

In many cases, the constant Ski,k2,k3,--X^) ™ ^q. ^ can be expressed in terms of natural numbers by identifying 
the corresponding series in eq. ^ when N = oo. Where this is not possible, it may be obtained numerically by 
comparing the expansion of eq. (|4]) in 1/N, once that has been determined by the method of this section, with the 
result obtained from eq. ([I]), for large, integer N. 
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IV. CALCULATION FOR LARGE |iV| WHEN ki = 1 



Some modification to the way the harmonic sum is calculated, eq. ([T]), is required before the above approach, namely 
the expansion of eq. Q at large |7V|, can be applied when ki — 1, because in this case Si^k2,k3,...{^) grows as (powers 
of) InA^ as increases. If /c2 7^ Ij the replacement 



N m N N N N m 

"2) = XI XI /("I' "2) + X /("I' "1) ~ X X /("2, "1) 

m — 1 n2 — 1 ni— 1 712 — 1 ni — 1 ni — 1 n2 — 1 



(15) 



in the definition in eq. ([3]) gives an expression from which a large |7V| expansion can be obtained: 



N 



ni — 1 



The expansion of the last term in eq. |T6|) is achieved by performing the usual analytic continuation of section |TT] and 
expansion of section Hill However, if fc2 = 1, then the replacement of sums in eq. (|15p can be used again in the last 
term in eq. (|16p . Repeating this process p times for a harmonic sum in which the first p indices equal 1 and fcp+i ^ 1 
gives 



p-i 



(17) 



where we define (we define the harmonic sum with no indices to be equal to one, i.e. S = 1) 



N 



KXJN) = (-1)" E ^-^^^^s^_Ms,,,,,,...M 

ni = l '^1 



(18) 



^kl k2 (-^) '^^'^ analytically continued and expanded at large \N\ according to the method of section [III] only 
when r 7^ 1, which excludes the last term in eq. (|17|) . Fortunately, we can derive from it the results we need to do 
this. Firstly, by performing the replacement of sums in eq. (|15p on the last term, eq. (|17p becomes 



p-i 



5w,fc..u...w=E(-i)""'^wW^w.fc.+.....w + E^r^^ 

I pP-l.sgn(fcp+i)(|fcp+i| + l)('JY\ I pP,kp+l 



(19) 



which can be analytically continued and calculated as a large \N\ expansion according to the method of section Hill 
This requires that the >5'i,....i(A^) for n — 1, . . . ,p are known in the form in cq. (These types of harmonic sums 

can also be determined by writing them in terms of harmonic sums of the same type but of lower depth according 
to formulae presented in Ref. 2]. This gives polynomial representations for 5*1,....! (-^) SiS determinants. However, for 

completeness we show how they can be determined within our approach.) To find these for n even, we use the result 
that follows from eq. ([T7]): 



n— 1 

Si,...AN) = ^ E 



771 — 1 



for n even. 



(20) 



5 



To obtain S'i,...,i(-/V) with n odd, we write it as X]ni=i - ^^lere S'i_...^i(ni) is expanded as in eq. (I20p 

with N ~+ m. Equation (|15|) is used repetitively on the second term, which then becomes 




(21) 



and eq. (fT5|) is used once on the first term in eq. ((20|) . The resuh is 




for n > 3 and odd. 



The quantity i^^'*^ (N) in eq. (|18p reduces to a harmonic sum when n — and/or the indices fci, fc2, . . . do not 
exist. Consequently, provided p is not too large, eq. (fT9|l can be used to easily express S'i,...,i,fcp+i,... (iV) in terms of 

_ p 

other harmonic sums as is done in Ref. |2f . For example, for p = 2, omitting the argument N for brevity, 

Sl.l,k = SiSi,k — <5'l,lS'fc + 5*2, A; — 'S'sgn(fe)(|fc|+l),l + ^kS,!- (23) 

For large values of p, this expansion into other harmonic sums becomes more difficult. In any case, this is not the 
goal of this letter. 



CALCULATION WHEN liVI IS NOT LARGE 



If jA^I is not sufficiently large for the evaluation of a harmonic sum as a large |7V| expansion to be sufficiently 
accurate, the harmonic sum can be calculated in terms of harmonic sums whose arguments are large. We begin with 
the relation between Sk-i.k2.k3,--X^) Sk-i,k2.ka.---i-^ + where r is any integer chosen such that |A'^ + r\ is large, 
that follows from eq. ([3]): 



Ski,k2M.---i^) = Sk,,k2M,--XN + r) - 



(sgn(fci))"^+^ ^ 

^'3/C2,/C3,...('^l +^)) 



1 — 1 ^ ^ 

Applying eq. (fM)l to Sk2M,--X'^i + then to 5'a;3,...(^i + and so on gives 



Sk,....,kAN) = Sk,....,k^N + r)- J2 



(sgn(fci))"i 



ni+N r 



'■-"1 f /I. \\ni+n2 + N r 



E 



(sgn(fc2))^ 



^^(ni+ 712 + AT) 1 



Sk2,...,kAN + r) 
Sks....,kAN + r) 



"2 = 



-N 



{ni + ... + np + 7V)l'=pl 



(24) 



(25) 
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The harmonic sums Ski,ki+i....{-^ + appearing in eq. ((25)) can all be calculated according to the method of section 
IIIII or lTVl i.e. as an expansion in 1/{N + r). 



VI. SUMMARY 



Our approach for evaluating a harmonic sum for any complex value of its argument can be summarized as follows. 
We first consider the case of large |iV|. When ki ^ 1, S'/ci,fc2,fc3,...(-^) ^^le form in eq. ^ can be obtained from 
Sk2,k3,...{^) the form in eq. ^ by analytically continuing Ski,k2,k3,...{^) to complex N using eq. (g]), writing 
Sk2M,---i^i + ^) ill the form in eq. ^ with N ^ ni + N , making the identifications in eq. (|lip . and finally using the 
expansion in eq. (fT2ll to put Ski,k2,...{N) in the form in eq. ([2|). Since the Sk{N) are known in the form in eqs. ([8|) 
and (fTO| . where ip{N + 1) has the large |A^| expansion in eq. ([6]), this means that all harmonic sums for which fci ^ 1 
can be put in the form in eq. Harmonic sums for which ki = 1 can be put in the form in eq. ([2]) in a similar 
manner. Si,,,,^i^kp+i....{-^) can be written in the form in cq. (jl9p . where Fj^^^f.^ (N) can be analytically continued 
p 

and expanded in 1/iV in precisely the same way as Ski,k2,k3,...{-1^) for fci 7^ 1 was, and can therefore be written in the 
form in eq. This requires that the S'i,...,i(iV) for n — 1, . . . ,p are known in the form in eq. which are found 

using eqs. (|20|) and (f22|) . When |A^| is not large, 5fci,fc2,fc3,...(-/V) can be calculated from harmonic sums with arbitrarily 
large argument using eq. (|25p . 

For convenience, FORTRAN and Mathematica files for calculating all harmonic sums of weight < 5 to double 
precision accuracy are provided at http://www.desy.de/~simon/HarmonicSums, which are suitable for leading order 
(LO), next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) calculations in McUin space. The 
extension of these programs to higher weights and higher accuracies is obvious. 

A preprint recently appeared [19i] which derives relations between harmonic sums of different weight. This allows 
harmonic sums up to a given weight to be calculated in terms of a few basic functions, each of which can be represented 
as a large |iV| expansion. 
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APPENDIX A: EXPLICIT RESULTS FOR LARGE |Ar| 

In this appendix we present explicit results for each harmonic sum with weight less than or equal to 3 as series in 
the inverse of its argument up to 0{1/N^'^) determined using the method in this letter. We define ln(7V) = ln(iV) + 7^ 
and 'r]{N) = (—1)^. The constants Ci_±i^±i appearing in Si^±i^±i{N) can be evaulated numerically as described in the 
last paragraph of section IIIII For example, they are given to 5 decimal places by ci_„i.i ~ —0.30883, ci^i._i — 0.51459 
and Ci__i._i = —0.37070. The constants S'-i,±i,±i(cx)) may also be evaluated numerically this way, and one finds 
S"-! _i'i(oo) ^ 1.47800, S"-!,! _i(cx}) = 0.66484 and S'_i,i,i(cx)) = -0.53721. Both the constants ci,±i,±i and 
S-i ±1 ±i(oo) may also be expressed in terms of rational and special numbers by using the results obtained in Ref. 

i-' ' 

1. Weight 1 



11 1 1 17 31 
S-i(N)^n{ ^ + r ^ + ^-^i-_)-ln2 (Al) 

Si{N) = InN + ^ - + - + - ^^^^ 



2. Weight 2 



TT^ 1 1 . . 

S2{N) — — — + ~ gjyS ~^ on ATF, ~ Am\T7 + onArfl (A4) 



1 


1 


1 


~^ 307V5 


427V7 


^ 30iV9 


31 ■ 


\ ln2 2 




' 4Ario. 


2 


+ 12 


1 


23 


9 


327V6 


168iV7 ^ 


1287V8 



2N 8N^ AN^ 16iV8 'liv-/ z iz /^^n 

1 3 5 1 19 1 23 9 263 19 



+ 



5_i,i(7V) ^r, 



2N 8N^ 24iV3 327V4 240iV5 32N^ leSN"^ 128N^ A80N^ 64iVio 

/I 7 1 61 73 379 11461 2041 371417 \ 

V27V2 " 24/V3 " 6iV4 + 240iV5 + 160A^6 ~ 504iV7 ^ 5040iV8 + isOiV^ + 20160iVio ) 

f 1 1 1 1 17 31 \ ,~ ^ ln2 2 7r2 



(A6) 



5i,_i(iV) 



2iV 4Af2 4^6 167V8 47Vio / 2 12 

/ 1 3 3 5 15 7 147 153 1185 



V4iV2 iQ]\j4 iQ]y5 327V6 8N'^ 647V8 32Ar9 64iVi 

~2]V + 12772 - + ~ + 1327V1" ' " 



(A7) 



Si i(N) = h 7T H r ^ H ? T7T lii^ H 



2 \2N 12N^ 120A^4 252iV6 2407V8 1327V1V 12 
1 3 1 1 1 1 1 221 3 23 



(A8) 



27V 8iV2 8iV3 288A/'4 48iV5 1440iV6 72iV7 604800iV8 iqqn^ 60480iVio 
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3. Weight 3 



fj. 3_ _5 2]_ 153 \ 3C(3) 

b-3[^) -Vy^^3 + 47V6 4Ar8 + 4Ario ) 4 



2m 2iV3 2N^ 2Af7 2iVV 4 
1 1 1 7 5 23 1 13 97 



4iV2 3iV3 47V4 i20Ar5 487V6 336iV7 goTV^ 160iVio 



11 151 469 331 67379 9181 



V2iV3 i2iv4 24iV5 240iV6 240iV7 126iVS 50407V9 4807Vio 

( _i l_ J. 3_ ~ 5C(3) 

+ ''V2iV2 2iV3 + 2Af5 2Ar7 + 2AfV « 

/ 1 _ 1 1 11 _ 13 _ 11 39 629 \ 

^ VTV 27V2^67V3 307V5 + 427V7 ) T^^~ 

52,i(Ar) = (-l + ^-^ + ^-^ + ^)rnAr + 2C(3) 

1 1 13 1 1 11 601 11 247 37 



N 4A/'2 367V3 gA/'* lOOA/'S 240iV6 352807V7 252iV8 75607V9 480iVio 

7rMii2 13C(3) / TT^ TT^ tt' tt' 177r^ 3l7r2 \ 

"^'"^^ ^ ~ 6 8 V 247V ^ 487V2 967V4 + 43^6 1927V8 + 487Vio ) 



1 7 5 19 2 173 

+ 



4Af2 12Ar3 8iV4 15Af5 24iV6 A2N^ 24N^ 5N^ 40N^° 

(m- f 12 + 3 22 - 7r2 1 77 - Stt^ 

-i,2( ) - 247V2 + 4/V3 " 48iV4 ~ 3^5 ^ 120^6 

31 1772 - 1197r2 206 2297 - 1557r2\ 7r2ln2_ 
30iV7 672iV8 35iV9 ^ 120iVio ) ^ ' 



/ 1 _ 5 1 7 _ 2 _ 27 27 187 \ _ Tr^lnTV C(3) 
Ai,-2( ) \^i]v^ - ^ + ^ + ^ - ~ 87V8 ^ 27V9 ^ 87^ 12 
n / 1 1 1 1 1 1 



24iV 144iV2 1440iV4 3024Ar6 2880iV8 1584Ario 



TT^lniV . 1 + -| - w . 17 , -150 + 77^ 



gi,2(iV)=^-C(3) + ^ + ^^ + ^+ ,,,,,, + 



TV 7V2 ' 367V3 720Ar4 450iV5 

7 TT^ QQ 3 I TT^ -1 -1 11 TT^ 

, 120 1512 _ I 56 1440 I _|_ 120 792 

2205m N» 350iV9 A^io 



(A9) 



Ss{N) = a3)-^ + ^-^ + ^-j^ + ^^ (AlO) 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 
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